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Abstract. The luminosity function (LF) of young star clusters is often approximated 
by a power law function. For clusters in a wide range of galactic environments this has 
resulted in fit indices near —2, but on average slightly steeper. A fundamental property 
of the —2 power law function is that the luminosity of the brightest object (£max) scales 
linearly with the total number of clusters, which is close to what is observed. This 
suggests that the formation of Young Massive Clusters (YMCs) is a result of the size of 
the sample, i.e. when the SFR is high it is statistically more likely to form YMCs, but 
no particular physical conditions are required. In this contribution we provide evidence 
that the LF of young clusters is not a —2 power law, but instead is curved, showing a 
systematic decrease of the (logarithmic) slope from roughly —1.8 at low luminosities to 
roughly —2.8 at high luminosities. The empirical LFs can be reproduced by model LFs 
using an underlying cluster IMF with a Schechter type truncation around w 2 x 
10^ Mq. This value of A/* can not be universal since YMCs well in excess of this Af^ 
are known in merging galaxies and merger remnants. Therefore, forming super massive 
clusters (> 10^ M0) probably requires conditions different from those in (quiescent) 
spiral galaxies and hence is not only the result of a size-of-sample effect. From the 
vertical offset a cluster formation efficiency of ~I0% is derived. We find indications 
for this efficiency to be higher when the SFR is higher. 



1 Introduction: the luminosity function of young star clusters 

It has been well established that the luminosity function (LF) of young star clusters is 
profoundly different from that of old globular clusters (GCs). Where the latter is peaked 
(when presented in logarithmic luminosity bins or in magnitude bins), the former rises 
down to the detection limit. When approximating the LF by a power law function of 
the form 
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dT 
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ma ny studies find a ^ 2 for clusters in wide ly varying environments (for overviews 



see lde Grijs et al.ll2003hlWhitmordl2003hlLars en 2006). H owever, the values foiind for 
a are often slightly higher (i.e the L F is steeper, see e.g. IDolphin & Kennicuttll2002l : 
Elmegreen et al.ll2002l : LarsenI 12002 ) and the range is rather large (L8 < a < 2.8) to 
be explained by statistical fluctuations or photometric uncertainties. 

The LF is often related to the cluster initial mass function (CIMF). However, the 
relation between the LF and the CIMF is not one-to-one, since the LF consists of clus- 
ters of different ages and evolutionary fading causes clusters of the same mass, but 
with different ages, to contribute at different luminosities in the LF. The empirically 
determined CIMF seems also well described by a power law, with an index of —2 (e.g. 



I 



lElmegreen & Efremovl[l997l : IZhang & Falilll999l : lMcCradv & Grahamll2007h . with the 
spread between different studies much smaller than that found for the LF. 

Although the LFs of the individual cluster systems are most of the time within 
a few a compatible with a —2 power law, a comparison between the different results 
shows that the deviations from the —2 power law are not randomly scattered, but instead 
show some systematic variations, in the sense that the LF is 



1 . steeper at higher luminosities ('Whitmore et al."l999';' 'Benedict et al.'2002l: lLarseril 
[2002 : .Mengel et al. 2005 : Gieles et al. 2006a ,b: Hwang & Lee 2008 ) ; 



2. steeper in redder fi lt ers (e.g. iDolphin & Kennic uttll2002l: lElmegreen et al.ll2002l : 



Sleeper in reuuer ji u ers le.j:^. | i^uiiJiim oc rs.eiimu uLuiz.uuzi. \cw 
Gieles et al.ll2006al lbl: iHaas et al.ll2008l : ICantiello et al.ll2009h . 



The steepening with increasing luminosity is illustrated in Fig. [J where we show 
results of power law fits to LFs of clusters in different galaxies taken from literature. 
The horizontal bars indicate the fit range used and the vertical bars give the uncertainties 
in the power law index. From this a decrease of —a (i.e. steepening) with increasing 
luminosity is seen. This suggests that the CIMF is not a universal continuous power 
law. If it was a continuous power law with the same index at all ages and in all galax- 
ies, all LFs should be power laws with the same index. Age dependent extinction or 
bursts in the formation rate would not cause a difference between the CIMF and the 
LF. An addition of identical power laws always results in the same power law. Lumi- 
nosity dependent extinction could cause a deviation from a —2 power law. This will be 
addressed in Section |3] 

The fact that the LF gets steeper at the bright-end is an indication that the CIMF 
is truncated at some upper mass A/up- When an abrupt truncation of the CIMF is 
assumed, it is possible to roughly estimate the index of the bright-end of the LF. Assume 
that the CIMF is fully populated, i.e. the mass of the most massive cluster actually 
formed, Mmax> is equal to A/up- Then assume a constant formation history of clusters: 
dA^up/dr =constant. The luminosities of these clusters, Lup, are age dependent. The 
light-to-mass ratio, or the flux of a cluster of constant mass, scales roughly with age 
as r~^^, with 0.6 ^ Ca ^ 1 depending on the wavelength A, such that dL^p/dr oc 
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. The same arguments hold for the luminosities of the 2"'^, 3'^, etc. 



most massive clusters, such that the bright-end of the LF is 
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For the F-band Cv = 0.7 (iBoutloukos & Lamersl2003l : lGieles et al.l2007l) . such that for 
the F-band LF the index at the bright-end is —2.5. For the U{I)-hmd Cu{i) — 1(0.6) 
and then equation ^ predicts a logarithmic slope of — 2(— 2.7), i.e. shallower(steeper) 
at blue(red) wavelengths. The faint-end of the LF should still be a power law with index 
—2 at all wavelengths or flatter if mass depe ndent disrupt ion is important. This double 
power law shape for the LF was found by iGieles et all f2006h) and , Whitmore et all 
(1999) for the LFs of clusters in M51 and the Antennae, respectively. 

In Section |2] we will model the LF in more detail and show that also this double 
power law shape is not the best approximation of the LF. 
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Figure 1. Sample of published indices of power law fit results to LFs o f youn 
star cl usters as a function of the fit range. The results are taken from: 'Larse_ 
(I2002h: six spiral s and the LMC; [M ora et al. (200l : NGC 45;|Gi eles et al. (2006b) 
iHaas et al. (2008) and Hwang & Lee (2008): M5 1 : 1 Whitmore et alj ( 119991) : Anten 
nae: , Dolphin & Kennicutt(2002.) : NGC 3627. 



2 Modelling the luminosity function 

In order to understand the luminosity function (LF) of star clusters we create semi- 
analytic cluster population models. For this the distribution function of initial masses 
(Mi) with age needs to be defined first. We define the probability of forming a cluster 
with an initial mas s between Mj ari d + dMj at a time between t and t + dt, as 



dN/{dMidt). For a lSchechteil (1 19761) type CIMF this then gives 



j^^ = AAf.-»p(-Af./M.), (4) 

where ^ is a constant that scales with the cluster formation rate (CFR) and is the 
mass where the exponential drop occurs. Note that we fix the power law index in the 
Schechter CIMF to -2. 

An equivalent distribution function for the present day masses (M) of clusters, 
dA^/(dMdt), can be found from multiplying equation ^ by \dMi/dM\ which de- 
scribes the relation between the initial and th e present cluster mas s to t aken into ac- 
count cluster dissolution. This was done by iFall & ZhangI (1200 ll) and I Jordan et"aL 



(120071) to mo del the mas s function of globular clusters using a constant mass loss rate. 



Lamers et al. I (12005) and lGielesI (120091 ) apply a mass dependent mass loss rate to model 



age and mass distributions of cluster populations using a pow er law and a S chechter 
CIMF, respectively. We refer the aforementioned studies and iLarsenl ( 20091) for full 
descriptions of the formulae. 

An equivalent distribution function of luminosity (L) and t can be acquire d by 
multiplying dN/ (d Mdt) by the age dependent mass-to-light ratio (T) (see also iFalll 
120061 : lLarsenll2009l ). such that 
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Figure 2. To/? panel: Model for the LF based on a constant cluster formation rate, 
a Schechter cluster initial mass function with = 2 x 10^ and a moderate 
cluster disruption time (full line). The dotted line shows a power law with index —2 
for comparison. Bottom panel: the relative contribution of different age groups. The 
logarithmic age ranges are indicates. 



diV diV ^, , 

T(t). (5) 



dLdt dMdt 

Since T(t) needs to be taken from SSP models (we adopt the lBruzual & Chariot! ( 2003b 



models for solar metallicity and a Salpeter stellar IMF) from here on the results need to 
be solved numerically. 

The luminosity function (LF) is then found by integrating AN/ (dLdt) over all ages 

diV f AN ^ 

-At. (6) 
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The logarithmic slope at each L, which is a probe of the shape of the LF, is found from 

dlog(djV/dL) 
dlog L 

An example of a model LF based on a Schechter CIMF and mass dependent dis- 
solution is shown in the top panel of Fig. |2l The bottom panel shows how clusters of 
different ages contribute to the LF. A feature of the truncated CIMF is that at the bright- 
end mainly young clusters contribute to the LF. This feature was indeed found from age 
determinations of the brightest clusters in a sample of spiral galaxies (Larsen,.2009) . 



3 Comparison to the cluster populations of M51, M74 and MlOl 

Here the model of the previous section is compared to the LFs of star clusters in 3 spiral 
gala xies: M5I (NGC 5 1 94), M74 (NGC 628) and M lOl (NGC 5457). We use the data 
from lHaas et al.l (l2008h : lH"wang & Led (l2008h (MSD. lBarmbv et al.l(l2006h (MlOl) and 
Gieles et al. (2010, in prep) (M74). 
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Figure 3. Logarithmic slope (—a) vs. luminosity for the V^-band LFs of clusters 
in M51, M74 and MlOl (left) and for 3 different filters for the M51 LP {right). The 
logarithmic slopes from the model shown in Fig.|2](using equation|7| are overplotted. 



First we take the model from Fig.|2]and apply a vertical offset such that it matches 
the empirical LFs. This directly gives us the cluster formation rate (CFR) which enters 
in the model through the variable A in equation Q. A comparison to the global SFR 
in these galax ies (Table 1) shows that the cluster formation efficiency (T = CFR/SFR, 
Bastian is roughly 10%, with a tendency of F to increase with increasing SFR. 



Note that this fraction could be affected by disruption if clusters are dissolving but are 
still picked up and therefore part of the LF. 



Table 1 . Comparison of the CFR derived from the LF of Fig.|2]to the SFR for the 
cluster populations considered here. 



Galaxy 


SFR 


CFR 


r 




[Moyr-i] 


[Moyr-i] 




M74 


-1 


0.05 


0.05 


MlOl 


-3 


0.25 


0.08 


M51 


-5 


0.90 


0.18 



Then we break the LFs in different luminosity segments and determine the loga- 
rithmic slopes in these intervals using a maximum likelihood estimate. The results are 
shown in Fig. |3l As can be seen the steepening of the LF with increasing L as shown 
in Fig. [T] is also found within the LF of individual cluster populations (left panel). In 
the right panel we show that the filter dependent logarithmic slopes for the M51 LF are 
also in good agreement with the model. 

So far we have ignored the effect of extinction on the individual clusters. If lumi- 
nous objects are more extincted, then the LF should be steeper in the bluer filter. Since 
the opposite is found, it seems that extinction does not play a vital role in shaping the 
LF 
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4 Summary 

We show that the LF of young star clusters is not a universal —2 power law, but instead 
shows evidence for curvature, in the sense that the logarithmic slope (—a) decreases 
towards brighter luminosities. This observational trend can be reproduced by a model 
based on a Schechter type cluster initial mass function with 2 X 10^ M©. A 

property of this model is that —a{L) is smaller (steeper) in redder filters at the same L, 
which is what we find for the rich cluster population of the interacting galaxy M51. 

Acknowledgments. MG thanks Pauline Barmby for kindly providing the lumi- 
nosities of the M 101 clusters and the organisers of the meeting Galaxy Wars: Stellar 
Populations and Star Formation in Interacting Galaxies for an enjoyable conference. 
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